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Multiply nested finite harmonic sums S ai ...a„(N) occur in many single scale higher 
order calculations in Quantum Field Theory. We discuss their algebraic and structural 
relations to weight w=6. As an example, we consider the application of these relations 
to the soft and virtual corrections for Bhabha-scattering to 0(a 2 ). 

1 Introduction 

Single scale processes in massless Quantum Field Theories [1,2], or field theories being 
considered in the limit m 2 /Q 2 — » [3], both for space- and time-like processes, exhibit 
significant simplifications when calculated in Mcllin space if compared to representations 
in momentum-fraction x— space. Here, the Wilson coefficients and splitting functions are 
expressed by Nielsen-type integrals 

f_1 )n+p+q-l rl J 

S n , P , q (x) = [ , / -ln( n - 1 )(z)ln^(l-^)ln«(l+^) 

1 [n)p.q. Jq z 

or harmonic polylogarithms [4]. The simplification is, to some extent, due to the structure of 
Feynman parameter integrals which possess a Mellin symmetry. The respective expressions 
can be expressed by finite harmonic sums S ai ... a „(N) for processes to 3-loop order [5,6], 
which form the appropriate language. Within the light-cone expansion, or analogous for- 
malisms for time-like processes, the argument of the harmonic sums are even- or odd inte- 
gers, depending on the process. However, one may consider mathematical generalizations, 
continuing the argument analytically to rational, real and complex values N e Q,R, C, re- 
spectively, [5,7]. In these extensions new relations between the harmonic sums are obtained, 
which lead to more compact representations. Since the hard-scattering cross sections usu- 
ally have to be convoluted with non-perturbative parton densities, it is convenient to widely 
work in Mellin space using analytic representations, also for the solution of the evolution 
equations. This also applies to the treatment of heavy flavor contributions in the full phase 
space, for which concise semi-analytic representations were derived [8]. The final x— space 
results are obtained by a single numerical Mellin-invcrsion performed by a contour integral 
around the singularities of the problem. 

In this note we give a summary on the algebraic and structural relations for finite har- 
monic sums, occurring in hard scattering processes. As an example we consider the virtual 
and soft QED corrections to Bhabha-scattering to 0(a 2 ) in the on-mass-shell scheme [9] to 
show that also this process fits to the general basis-representation being derived for various 
other two- and three-loop QCD processes. 
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2 Algebraic Relations 



The complexity of finite harmonic sums is given by N max = 2 ■ 3 W , with the weight 
w = X)fc=i l°fcl growing exponentially In the x— space representation of the 2-loop Wilson 
coefficients in QCD, which are described by w = 4, nearly 80 functions emerge, which 
corresponds to the maximum level possible, cf. [5, 10]. The algebraic relations of finite 
harmonic sums operate on their index set and are implied by their quasi-shuffle algebra [11]. 
The algebraic relations of finite harmonic sums were investigated in [12] in detail. An 
example for the shuffle product is 

Sa± ,a,2 L- fJ ^ 0,3,0,4 — S ai , a2 , a3 , a4 -\~S ai ,0,3,0,2 ,04 ~\~S ai ; a 2 jft4 ; a 2 ~\~S a3j a 4 , ai , a2 -\~S a3 ) a li fl 4; tt 2 ~\~S a3 , ai ;fl2 ;tl4 . 

For a given index set the number of basic harmonic sums is counted by the number of Lyndon 
words and can be calculated by a Witt-formula. Investigating the type of harmonic sums 
emerging in physical single scale problems up to 3-loop order, cf. [1-3], the index { — 1} 
never occurs. Their number is iV-,_i = [(1 — \/2) w + (1 + \/2) w ]/2 [13]. One therefore may 
significantly reduce the number of basic functions using N-,-i in the corresponding Witt 
formula. The following table illustrates the corresponding complexities of combinations # c 
and the number of algebraic basis elements #5 in dependence of the weight w. The initial 
complexity in case of the absence of indices a; = — 1 is lower than the number of algebraic 
basis elements in the complete case. 



w 


1 


2 


3 


4 


5 


6 


#c 


2 


8 


26 


80 


242 


728 


#b 


2 


5 


13 


31 


79 


195 


#c(- - 1) 


1 


4 


11 


28 


69 


168 


#6(-"l) 


1 


3 


7 


14 


30 


60 



3 Structural Relations 



The algebraic relations lead to a first reduction of the complexity of finite harmonic sums. 
They are independent of the value of these quantities. Beyond these relations, the structural 
relations are of a more specific type, cf. [5,14]. There are three types of these relations. The 
first class is implied by allowing half-integer values for N. A second class emerges through 
partial integration using the representations of harmonic sums through Mcllin-transforms 
e.g. of weighted harmonic poly logarithms. The third set is implied by differentiating har- 
monic sums w.r.t. their argument, which requires JVeR. To illustrate case 1 we represent 

1 1 

1 -x 2 ~ 2 

The Mellin transform of this equation implies, that S^i(N) is linearly dependent of Si(N), 
if N e Q. Various relations of this type emerge at higher weight [14]. If N e R, one may 
differentiate harmonic sums and obtains 

^-S ai ... „(iV) = J' dxx N ~ 1 Hx)f(x) . 

S ai ...a„(N) is given by the Mellin transform of f(x). One may represent the derivative 
in terms of polynomials of harmonic sums and their values at N — > 00. The analytic 
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continuation of a harmonic sum to N e C is a meromorphic function with poles at the non- 
positive integers. Up to polynomials growing oc S™(N), \N\ — > oo, which can be separated, 
harmonic sums are factorial series. I.e. they obey an analytic recursion relation F{z + 1) — ► 
F{z) and one may calculate their asymptotic representation for z — ► oo analytically. In this 
way these functions are uniquely defined in the complex plane. Due to the above properties 
we now define equivalence classes of representations, which contain a harmonic sum and 
all its derivatives. We will only count these equivalence classes, since all derivatives can be 
easily found analytically if the respective lowest weight sum of the class is known. In the 
following table we summarize the functions of the basic Mellin transforms up to w = 5. The 
case w = 6 is still to be completed. 



w = 1 l/(x - 1)+ 
w = 2 ]n(l + x)/(x + l) 
w = 3 Li 2 (x)/(x±l) 
w = 4 Li 3 (V)/(a; + l) 
w = 5 U 4 (x)/(x±l) 

Li§(s)/Cr + 1) 
w = 6 Li 5 (x)/(.T + 1) 

S 3 ,2(a:)/(a:±l) 



S h2 (x)/(x±l) 

S h3 (x)/(x + l) S 2 , 2 (x)/(x±l) 

[S 2 ,2{-x)-Ul{-x)/2]/{x±l) 

S 1A (x)/(x±l) S 2 , 3 (x)/(x±l) 



These functions emerge in the following processes : number of fct. 



• 


0(a) 


w 


< 2 


Wilson Coefficients/anom. dim. 


#1 


• 


0(a 2 ) 


w 


< 3 


Anomalous Dimensions 


#2 


• 


0(a 2 ) 


w 


< 4 


Wilson Coefficients 


# <5 


• 


0{a 3 ) 


w 


< 5 


Anomalous Dimensions 


#15 


• 


0{a 3 ) 


w 


< 6 


Wilson Coefficients 


#29+ 



4 Example: Bhabha Scattering 

As an example we consider the soft- and virtual corrections to Bhabha-scattering up to 
O(o 2 ), T 2 (x) = S^\x) x (1 — x + x 2 )/x 2 , [9]. The x— space expression simplifies in Mellin 
N space and depends on three basic functions £2,1,1, £2,1 and Si only aside of polynomial 
pre-factors : 



248 -|- 15 N^ 1 -f- 7V^ 2 
T2{N) = 2(JV - 2)(JV - 1)JV(JV + 1)(N + 2) Sl ^ N) ~ (JV - l)(JV + 1) S *^W 

-340 + 120 AT + 17 N 2 + 18 JV 3 - 31 Af 4 , N -304 + 278AT + 81AT 2 - 38Af 3 + 19AT 4 , , 

S3,l(ff) — — — -— — ; — — ; oX Sa{N) 



2(N -2)(AT -1)N(N + l)(A r + 2) ' 8(AT - 2)(N - 1)N(N + 1)(N + 2) 

304 - 328 N - 500 N 2 + 330 N 3 — 6 + 6 — 2 + 4 N 7 , 

S2 1 (N) 

(N-2) 2 (N-l) 2 N 2 (N + l)(N + 2) ' K ' 



(N -2)(JV -1)N(N + 1)(N + 2) ' v ' v ' (JV — 1)N(N + l)(A r + 2) 
-1840 + 292 N + 5532 TV 2 + 827 TV 3 - 1978 TV 4 - 274 TV 5 + 36 N 6 + 19 N 7 - 22 N s 
i(N - 2) 2 (N - 1) 2 N 2 {N + 1) 2 (N + 2) 



S u ,i(W) 
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128 - 56 N - 252 A" 2 + 54 N 3 + 177 A" 4 — 91 TV 5 + 19 A" 6 + 9 TV 7 

2(iV - 2)(AT - 1) 2 N 2 (N + 1) 2 (N + 2) 3( ' 

4032 - 2048 N - 14200 iV 2 + 5036 AT 3 + 23610 N 4 + 2521 - 12342 N 6 

4(N - 2) 3 (N - 1) 3 N 3 (N + 1) 3 (N + 2) 
-3365 N 7 + 2148 N 8 + 903 TV 9 + 14 A" 10 - 167 AT 11 + 50 N 12 
4(N - 2) 3 (N - 1) 3 N 3 (N + 1) 3 (N + 2) 



Si,i(N) 



Si,i( N ) 



124 + 16 N + 24 N 2 - 4 N 3 - 14 N 4 , s , , 424 — 118 A" + 9 A" 2 — 2 A 1 " 3 + 23 N 4 , , , v 

Si,i(iV)C(2) + / —4 oN S 2 (N)Si,i(N) 



(N - 2)(N - 1)N(N + 1)(JV + 2) ,v / w 4(AT — 2)(A r — 1)N(N + l)(A r + 2) 

224 + 144 Af - 1216 N 2 - 56 AT 3 + 1786 N 4 + 641 A" 5 - 406 N 6 

4(7V - 2) 2 (AT - 1) 3 AT 3 (AT + 1) 3 (7V + 2) 2( ' 

17 N 7 - 308 N s + 141 AT 9 - 56 AT 10 + AT 11 58 + 21 N + N 2 + 15 A^ 3 + 10 A^ 4 rjvv ,„, 

4(Ar-2) 2 (AT-l) 3 A r3 (Ar + l) 3 (Ar + 2) 2( (AT — 2)(AT - l)JV(iV + 1)(AT + 2) 2( )C( ' 

232 - 384 N 2 - 17 N 3 + 286 N 4 - 128 A" 5 — 14 A" 6 + A" 7 

4(AT - 2){N - 1) 2 N 2 (N + l) 2 (Af + 2) 2( ) i( ) 

-560 — 26 A" — 31 A" 2 — 10 A^ 3 — 33 A" 4 . >tn2 

#2 (AO 

8(AT - 2)(JV - 1)AT(7V + 1)(AT + 2) v ; 

576 + 1088 N - 3280 N 2 - 5136 N 3 + 11764 N 4 + 20392 N 5 - 17385 AT 6 - 30114 N 7 

4(AT - 2) 3 (Ar - 1) 4 N 4 (N + 1) 4 (N + 2) 

5984 N 8 + 17228 A" 9 - 1228 N 10 - 2754 AT 11 - 112 A^ 2 - 8 TV 13 + 33 N 14 - 24 N 15 



Si(N) 



4(Af - 2) 3 (N - l) 4 N 4 (N + 1) 4 (AT + 2) 
-56 + 336 N + 522 N 2 + 424 N 3 - 53 N 4 - 500 AT 5 + 60 A" 6 + 28 A" 7 — 5 N 8 
2(AT - 2) 2 (iV - l) 2 Af 2 (Af + l) 2 (Af + 2) 



Si(JV) 



Si(JV)C(2) 



64 + 6 N 2 + N 3 2112 + 608 AT + 76 AT 2 - 140 N 3 + 107 N 4 N , 
Sl(NK(3) + -TTTr 77777 7777777 , T7777 , ^ C( 2 ) 



(N-2)(N- 1)N(N + 1) 10(JV-2)(JV- 1)JV(JV + 1)(JV + 2) 

-224 - 136 AT + 1688 AT 2 + 1290 N 3 - 1998 AT 4 - 1997 + 198 TV 6 

2(7V - 2) 2 (Af - 1)3N 3 {N + l) 3 (Af + 2) ^ (2) 

405 Af 7 + 376 AT 8 - 119 N 9 + 56 AT 10 + 5 A" 11 , , 

' C(2) 

2{N - 2) 2 (N - 1) 3 N 3 (N + l) 3 (Af + 2) S W 

-552 + 144 Af + 1654 Af 2 - 370 A" 3 - 361 N 4 + 19 A" 5 + 35 - 25 Af 7 



2(Af - 2) 2 (Af - l) 2 Af 2 (AT + l) 2 
Pi (AT) 



C(3) 



16(AT - 2) 3 (Af - 1) 5 (AT + 1) 5 N 5 (N + 2) 

N 4 - N 2 + 12 Af 4 - Af 2 + 12 



" T "* {N - 2)(N - 1)N(N + 1)(N + 2)^°' 2 2 (TV - 2)(JV - 1)N(N + l)(N + 2) 
with 

Pi = 320 - 64 N - 1920 Af 2 + 1600 Af 3 +6524 N 4 - 14872 N 5 - 19036 N G +31543 Af 7 
-43960 Af* - 13935 Af 9 + 65372 Af 10 + 26822 Af 11 - 44576 Af 12 - 9558 Af 13 
+9840 Af 14 + 339 Af 15 + 428 N 1(i - 371 N 17 + 128 Af 18 

There are no alternating sums contributing, unlike the case for other 2-loop Wilson coeffi- 
cients in QCD [1,3], where generally up to 6 functions emerge. 



5 Conclusions 

The single-scale quantities in Quantum Field Theories to 3 Loop Order, corresponding to 
w = 6, can be represented in a polynomial ring spanned by a few Mellin transforms of 
the above basic functions, which are the same for all known processes, including QED 
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processes, as the soft- and virtual corrections to Bhabha-scattering. This points to their 
general nature. The basic Mellin transforms are meromorphic functions with single poles 
at the non-positive integers. The total amount of harmonic sums reduces due to algebraic 
relations as a consequence of the index structure, and structural relations continuing the 
argument to N e Q, N e R. They can be represented in terms of factorial series up to simple 
"soft components" . This allows an exact analytic continuation. Up to w = 6 physical 
(pseudo-) observables are free of harmonic sums with index ={ — 1}. To w = 5 all numerator 
functions are Nielsen integrals. 
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